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The effective Hamiltonian for perovskite titanates is derived by taking into account the three-fold 
degeneracy of tig orbitals and the strong electron-electron interactions. The magnetic and orbital 
ordered phases are studied in the mean-field approximation applied to the effective Hamiltonian. 
A large degeneracy of the orbital states in the ferromagnetic phase is found in contrast to the case 
of the doubly degenerate e g orbitals. Lifting of this orbital degeneracy due to lattice distortions 
and spin-orbit coupling is examined. A general form for the scattering cross section of the resonant 
x-ray scattering is derived and is applied to the recent experimental results in YTiC>3. The spin 
wave dispersion relation in the orbital ordered YTi03 is also studied. 



PACS numbers: 75.10.-b, 71.10.-w, 78.70.-g, 71.30. 
I. INTRODCUTION 

Since the discovery of layered superconducting 
cuprates I^-xSr^CuO^ studies of electronic structures 
of transition-metal oxides are revived pfrom the mod- 
ern view point of electron correlation.^ Perovskite ti- 
tanates Ri- x A x TiOz are one of the prototypical three- 
dimensional materials which show the Mott transition 
and anomalous metallic states at a vicinity of the tran- 
sition. Here, R and A indicate the trivalent and divalent 
cations, respectively. The end compounds i?TiC>3, where 
a nominal valence of all Ti ions is 3+, are recognized to 
be Mott insulators. A mismatch of the ionic radius in the 
pseudo-cubic structure induces the GdFe03-type lattice 
distortion, i.e. a tilting of a TiOg octahedron.0 The elec- 
tronic structure of i?TiC>3 systematically changes with a 
kind of the rare-earth ion i?pa'Q a large distorted YTiC>3 
shows a ferromagnetic ordering at 29K. The saturated 
magnetic moment is 0.84jj which is close to the expected 
value from S — lpZ.2. and a definite optical gap is observed 
to be about 2eV.Lfl3 On the contrary, the insulating char- 
acter of the less distorted LaTiC>3 is more marginal than 
YTiC>3; a staggered magnetic moment for the G-type an- 
tiferromagnetic_(AF) state is less than one half of the 
expected valueo and an insulating gap is smaller than 
0.5eVJ3'Ei By doping i?TiC>3 with holes, the system un- 
dergoes the metal-insulator transition and exhibits sev- 
eral unconventionalj-picjtallic characters ascribed to the 
electron correlation jailil as well as La2- x Sr x Cu04. 

One of the remarkable discrepancies in perovskite ti- 
tanates from layered superconducting cuprates is that 
the orbital degree of freedom survives in titanates; the 
electron configuration of Ti 3+ is d 1 where the three t^g 
orbitals, i.e. d xy , d yz and d zx orbitals, are degenerate 
under the cubic crystalline field. Thus, this ion has a 
degree of freedom which indicates an occupied orbital by 
an electron. The intensive and extensive studies of the 
orbital degree of freedom have been carried out recently 
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in colossal magnetoresistive (CMR) manganitesEl where 
a Mn 3+ ion has doubly degenerate e g orbitals. Here, 
it is widely believed that observed dramatic phenomena 
such as CMR are caused by strong interplay between 
spin, charge and orbital as well as lattice. In comparison 
with manganites, there exist following characteristics in 
the orbital degree of freedom in titanates: (1) there is a 
three-fold degeneracy of the t 2g orbitals, (2) an electron 
hopping between nearest-neighboring (NN) different or- 
bitals is prohibited in a cubic crystal structure, (3) the 
spin-orbit (LS) coupling is possible to be relevant- and 
(4) the cooperative Jahn- Teller (JT) effect is weak.13 Ac- 
tually, some theoretical and experimental studies of per- 
ovskite titanates have been dpaeJjrom the view point of 
the orbital degree of freedom.E^TcS In addition, the reso- 
nant x-ray spattering (RXS) was applied to i?Ti03 very 
recentljiiSEj and the orbital ordering was successfully ob- 
served in YTi03t3 The systematic studies by utilizing 
this experimental method are expected to clarify roles 
of the orbital degree of freedom in the Mott transition 
and the several unconventional phenomena observed in 
titanates. 

In this paper, the effective Hamiltonian for the elec- 
tronic structures in perovskite titanates is derived and 
the spin and orbital structures in i?Ti03 are studied 
in the mean field approximation. The derived Hamil- 
tonian corresponds to the tj model in superconduct- 
ing cuprate£3 and the spin-orbital model for CMR 
manganites £j and is applicable to a wide range of doped 
and undoped titanates. In particular, we focus on roles of 
the three-fold degeneracy in the t^g orbitals. It is shown 
that the orbital degeneracy of the ground state is more 
significant than that in the e g orbital case. A general 
form for the scattering cross section of RXS is derived and, 
is applied to the recent experimental results in YTiOs.EHl 
Roles of the orbital on the spin wave dispersion relation 
in YTi03 are also discussed. 
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II. EFFECTIVE HAMILTONIAN 

We start with the tight-binding Hamiltonian in a three 
dimensional lattice consisting of Ti ions. Three ti g or- 
bitals d xy , d yz and d zx and the intra-atomic Coulomb 
interactions are considered in each Ti ion. The Hamilto- 
nian is 

H = £ 

(ij)-ry'o- 

+ Uy]ni^fm^i + U'- n n n n > 

27 

+ I d\ la d\^, a ,di 1<y i diy a 

*7>7' ' crcr' 

+ I £ 4^ d ni d ^'l^Yh 

«7#7' 



(1) 



where rft creates a ti g electron at site i with spin 
cr(=|,|) and orbital 7(= xy,yz, zx). U and {/' are the 
intra-orbital and inter-orbital Coulomb interactions, re- 
spectively, and / is the exchange interaction. In an iso- 
lated ion, these interactions are represented by the Racah 
parameters as U = A + 4B + 3C, U' = A - 2B + C and 
I = 3£> + C, and a relation U = U' + 21 is satisfied. 
tJJ is the hopping integral between site i and its NN 
site j with orbitals 7 and 7', respectively. In a simple 
cubic lattice, the hopping integral is diagonal and one 
of the diagonal components is zero. For a Ti-Ti bond 
in a direction l(= x,y, z), two orbitals, which have a fi- 
nite hopping integral tj^ 1 , are termed active orbitals de- 
noted by ai and and one with no hopping integral is 
termed an inactive orbital denoted by cj. For example, 
(a x , b x , c x ) = (zx, xy, yz). tj? is simply expressed in this 
case as 



77 _ 



t8 11 '{8 iai + 5 lbl ). 



(2) 



The GeFe03-type lattice distortion breaks this relation 
as discussed in Sec. [II. 



Since the Coulomb interactions U aiuL.U' are larger 
than the hopping integral in titanatespil the effective 
Hamiltonian is derived by perturbational calculation 
with respect to the hopping integral. The Hamiltonian 
is 



T~t — TLt + T~tj, 



(3) 



where the first and second terms correspond to the so- 
called t and J terms in the t J model, respectively. The 
t term is given by 



T\ 



B.c. 



(4) 



where d ila (= di 7CT II( 7 v')^( 7 ff) (1 - d] 1 , a ,d ll ' a ')) excludes 
multi occupied states of electrons at site i. The J term 



is classified by the point symmetry of the intermediate 
electronic states, i.e. d 2 states. In the case where Eq. (^) 
is satisfied, the Hamiltonian is 



with 
Ht 2 



'Hj = + 'Hti + He + 7~Cai , 



-^£ 

(ij) 



(5) 



Si'S j ][B l -C l + D l ), (6) 



- J T,Y.{\ n ^-^-^J ){B l + C l +D l ), (7) 



(ij) 



He = 



lifij Si ■ Sj 



-A 1 - -C n 
3 3 



(8) 



H * = - J ^ £ (i n w - $ ■ (l Al + l c ") ■ ( 9 ) 

Prefactors are given by J Tl = t 2 /(U'-I), J Ta = t 2 /(U' + 
I), J E = t 2 /(U - I) and J Al = t 2 /(U + 21). By using 
the relation U = U' + 21, we obtain Jt 2 — Je- n>i(= 
J2a-y d\ ia di ia ) is the number operator and Si is the spin 
operator given by 



2 £ 4 7 <7^' d n^'- 

■yaa' 



(10) 



A 1 , B l , C l , C n and D l are the orbital parts of the 
Hamiltonian represented by the eight orbital operators 
Or 7 where T denotes an irreducible representation in the 
Oh group and 7 classifies the bases of the irreducible 
representation. To represent Or~/, let us introduce the 
Gell-Maan matrices which are generators of the SU(3) 
algebra:E£l 
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We define the orbital operators as 



-1 



CJQ/3 



(15) 



where (T T ,l) ={Eu;S), (Ev,3), (T 2 x;6), {T 2 y;4), 
(T 2 z;l), (Tix;7), (Tiy;5) and (T x z;2). The operators 
OiE-j and OiT 2l describe the electric quadrupole mo- 
ments and OiTi7 describes the magnetic dipole moment. 
By utilizing the above operators, we obtain 
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C l — 2{OiT 2 lOjT 2 l + OiTilOjTil), 



(16) 
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FIG. 1. Schematic pictures of representative virtual ex- 
change processes for the terms (a) A 1 , (b) B , (c) C ! and (d) 
D l defined in Eqs. @, (0), @ and @, respectively, a; 
and bi are the active orbitals and c; is the inactive orbital (see 
text). 
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where 0' B7 's are given by 
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COS =£-771; Sin =£-777; 

-sin^rn/ cos ymj 
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(19) 



(20) 



(21) 



with mi = (1,2,3) for I — (x,y,z). It is worth rewrit- 
ing the orbital parts of the Hamiltonian from the view 
point of the active and inactive orbitals; we introduce 
the pseudo-spin operator with the quantum number 1/2 
for a bond along a direction I 



rpl 

7,7' = (a(,6i),cr 



1 



E 



(22) 



with c? being the Pauli matrices, and the number operator 

n -ni = Z)<7 ^ff^-yiff for 7i = (o/j &z, c/). The orbital parts 
of the Hamiltonian are rewritten as 

^' = 4(n iai nj ai +n. lbl n jbl ), 
B l = 2(n lai n jbl + n ibl n jai ), 

A ( rpl rpl | rpl rpl \ 

^ \ ix jx ' iy jy ) 3 

' / y| / rpl rpl rpl rpl \ 

^ ^\ J -ix 1 jx iy jy)' 
D l = n ici (n jai + rijb,) + (n iai + n lbl )n JCl . 



(23) 
(24) 
(25) 
(26) 
(27) 



We note that T['s with different I are not independent 
with each other. These simple expressions result from 
the diagonal form in the hopping integral under which an 
electron number of each orbital is conserved. Schematic 
pictures of representative exchange processes for A 1 , B l , 
C l and D l are shown in Fig. [I]. 

A similar model Hamiltonian with Eq. (|5|) was derived 
in Ref. ^?], although the intra-atomic exchange interac- 
tion I is assumed to be zero. As shown in the next sec- 
tion, this condition corresponds to a critical point in the 
phase diagram. Also, a similar model with / = is rep- 
resented by the pseudo-spin operators T- (Eq. (p2j) ) in 
Ref. [l9|. In a model Hamiltonian derived in Ref. go], two 
of the three t 2g orbitals are taken into account. 



III. SPIN AND ORBITAL STATES 

The effective Hamiltonian is applied to i?TiC>3 where a 
valence of all Ti ion is 3+. Consider a hypothetical-cubic 
lattice consisting of Ti ions, instead of the actual crystal 
lattice of RTiO^- The mean- held approximation at zero 
temperature is applied to the Hamiltonian and the four 
kinds of spin and orbital ordered states are considered; a 
uniform spin (orbital) state termed F and three staggered 
spin (orbital) states termed A-AF, C-AF and G-AF. The 
periodicities of these orderings are characterized by the 
momentum (000), (007r), (ttttO), and (tttttt), respectively. 
(Siz){= ±1/2) is adopted to be a spin order parameter 
and, the orbital order parameters (0^7) are calculated 
by the orbital wave function at site i 
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FIG. 2. (a): The spin and orbital phase diagram as a 
function of Rj(= Jax/Jtx)- Schematic pictures of the repre- 
sentative orbital ordered states for Rj < 1 and Rj > 1 are 
shown in (b) and (c), respectively. 

The order parameters are optimized numerically to ob- 
tain the lowest energy. In Fig. |^, the magnetic and or- 
bital phase diagram is presented as a function of Rj = 
Jai/Jt ± - The relation Jt^/Jt^ — 5Rj/(2 + 3Rj) derived 
from the condition U = U' + 21 is used. Although a 
value of Rj is smaller than one in actual compounds, the 
calculated results in the region of Rj > 1 are also shown 
for comparison. In the region of Rj < 1, the F spin 
state and the G-AF orbital ordered state are realized. 
The wave functions in the A and B orbital sublattices 
are given by 



1^ 



\d a ), 



\ip B ) = cos 6^) +sin(9e l<# X 



(29) 



respectively, with any values of 8 <E [0, ir] and (j> € [0, 2ir], 
and (a, (3, 7) = (xy,yz,zx), (yz,zx,xy), (zx,xy,yz). 
There is a large continuous degeneracy in the orbital 
state; any linear combinations of \dp) and |<i 7 ) are de- 
generate in the sublattice B. This is generalized to the 
states 



\A A ) 

\4>Sb) 



= \d a ), 



cos I 



JB 



\dp) 



sin f 



\d 



111 



(30) 



where %a (is) indicates the i(j)-th site in the A(B) or- 
bital sublattice. 8j B and <pj B at each site are taken inde- 
pendently, because the Hamiltonian includes the interac- 
tions between the NN sites. In order to understand this 
state in more detail, let us consider the hopping integral 
between the NN occupied orbitals defined by 



4>i 



<jj>77'<7 



tl7'dld 



H.c. 



(31) 



Tjj's are zero for all bonds in this orbital ordered 
state. Therefore, the exchange processes denoted by A 1 



(Eq. (16)) do not occur (see Fig. and the ferromag- 
netic interaction is dominant. It is worth comparing the 
present results with those in the system where the dou- 
bly degenerate e g orbitals, i.e. the d 3z 2_ r 2 and d x i_ y i or- 
bitals, exist. Here, the effective Hamiltonian correspond- 
ing to Eq. (H|) is expressed by the spin operator Si and 
the pseudo-spin operator for the orbitaLdegree of free- 
dom Ti with the quantum number l/2.t3n3 The orbital 
state in the ferromagnetic phase obtained by the mean- 
field approximation is the G-AF orbital state where the 
wave functions are given by 



\*I>a) 
\1>b) 



■(1) 
*>(i)l 



3z 2 -r< 



'3z 2 - 



(32) 



for any value of G [0, 27r].S This is also a staggered or- 
bital ordered state with a continuous degeneracy. How- 
ever, there exist the remarkable differences between the 
t2 g and e g cases: (1) The hopping integrals between the 
occupied orbitals Ty are finite in the e g case and de- 
pend on 9. Thus, the generalization of the orbital or- 
dered state, as seen in the t 2g case (from Eq. ( p9| ) to 
Eq. (pp|)), is impossible. That is, the orbital degeneracy 
in the ground state is more remarkable in the ti g case. 
This is because a number of the orbital degree of free- 
dom is larger in this case. (2) The orbital wave functions 
with complex coefficients are not included in the e g case. 
This is because the effective Hamiltonian for the e g elec- 
tron is represented by the operators T x and T z , unlike T y 
which breaks the time reversal symmetry. In the region 
of Rj > 1 in Fig. 0, there is no continuous orbital de- 
generacy. The C-AF and G-AF spin states are realized 
associated with the F orbital state with \tp) — \d a ) for 
a = (xy,yz,zx). The point Rj = 1, i.e. I — 0, is a 
critical point in the phase diagram where degeneracy of 
the spin and orbital states are significa|H±. It is supposed 
that Rj for i?Ti0 3 is about 0.2 - 0.4P although there 
is a large ambiguity in an estimation of I. 

We next examine how this large continuous orbital de- 
generacy is lifted by the following three effects observed 
in i?TiC>3: the GdFe03-type lattice distortion, the JT- 
type distortion in a TiOg octahedron and the LS cou- 
pling. The GdFeO-3-type lattice distortion bends a Ti- 
O-Ti bond. The simple form of the hopping integral 
in Eq. (^|) is not valid and tj^ for any pairs of 7 and 
7' are finite. We calculate all components of the hop- 
ping integral ij^ for the crystal structures in YTi03 
and LaTiOs by the Slater-Koster formula.cHl The most 
remarkable changes are found in the hopping integrals 
between the different active orbitals t""' b ' . This is be- 
cause the GdFe03-type distortion induces a a bond be- 
tween the d ai {db t ) orbital and the 2pi orbital at the NN 
O site. We simulate this distortion by introducing a new 
term in the hopping integral as 



t 



(33) 
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FIG. 3. A schematic picture of the crystal structure of 
iiTiOs. Filled and open circles indicate Ti and O ions, re- 
spectively, and Ai, A2, B\ and B2 denote Ti sites in a unit 
cell. Arrows indicate displacements of O ions in the ab plane, 
and crosses and dots in the O ions indicate the displacement 
along the +c and — c directions, respectively. 

where the sign of the tranfer integral s; = ±1 depends on 
the direction I and a parameter R t = t a /t^ is interpreted 
to be an increasing function of this distortion. A value 
of Rt for YTIO3 is estimated to be about 0.5 ~ 1. The 
effective Hamiltonian including the GdFe03-type distor- 
tion, the JT-type distortion and the LS coupling is given 
by 

H = Hj + Hj T + H L s, (34) 

where t term is neglected. The first term is the modified 
J term including the GdFe03-type distortion and its ex- 
plicit form is presented in the Appendix [a|. The second 
term is for the JT coupling 

Wjt = 9e ^2 O 

+ 9T 2 ^ O l T 2l QiT- 2ll (35) 
i*y=x,y ,z 

where Qr-y's are the normal coordinates of the oxygen 
displacements in a TiO@ octahedron with the symmetry 
L7, and gr's are the coupling constants. The last term 
is for the LS coupling denoted by 

Ti-LS = Als ^ SijOiT^- (36) 

i*y=x,y ,z 

The corrections for the effective Hamiltonian in Eq. j3^ ) 
are of the order of 0(g r t 2 /U 2 ) and 0(X LS t 2 /U 2 ). In the 
mean field calculation, an orthorhombic unit cell with 



1.5 



0.0 

1.0 



(a) 






F 


\ A-AF 


G-AF 
C-AF 


.(b) 


C-AF \ \ 


G-AF 


r 


C/G-AF.\ 




. (c) 


t 1 1 1 1 


G-AF . 


F 


C/G-AF 


XyL 



$00 0.25 0.50 0.75 1.00 1.25 



FIG. 4. The magnetic phase diagrams as functions of 
Rj and (a) a ratio between the hopping integrals R t = t^/t^ 
caused by the GdFeOs-type lattice distortion, (b) the JT-type 
lattice distortion gQ, and (c) LS coupling Xls- Other parame- 
ter values are chosen to be Xls ~ Rt = in (a), gQ — Xls = 
in (b), and gQ — R t = in (c). 

the Pbnm space group is adopted. Four Ti ions in a 
unit cell are termed Ai, A2, B\ and B2 (see Fig. ||) 
and the orbital states at these sites are considered in- 
dependently. As for the JT-type distortion, Qr 2 y an d 

-j-Qeu + ^-Qev{= QE3x 2 -r 2 ) are dominant at site A\ 
in YTiC>3. Thus, we assume in this calculation that 
\9n\QAiT2y = \gE\QA 1 E3x 2 - r i(= gQ) and other com- 
ponents are zero. The JT-type distortions at other sites 
are introduced by considering the Pbnm space group. 

The phase diagrams with including the GdFeC>3-type 
lattice distortion, the JT-type distortion and the LS cou- 
pling are presented in Figs. |i| (a), (b) and (c), respec- 
tively. In all phase diagrams, the large orbital degeneracy 
is lifted and a region of the ferromagnetic phase shrinks. 
This result implies that the high symmetry in the orbital 
space is favorable for the ferromagnetic spin ordering. 
The characteristics of the each phase diagram are sum- 
marized as follows; (1) By introducing the GdFe03-type 
lattice distortion, the orbital degeneracy is partially lifted 
and one of orbital ordered states in the ferromagnetic 
phase is given by \ip Al ) = 0.71|da, y ) + 0.71\d yz ), \ip Bl ) = 
0.38\d xy ) - Q.35\d yz ) + 0.85|d z:c ), \i>A 2 ) = 0.38\d xy ) - 
0.85\d yz )+0.35\d zx ) and \yj B . 2 ) = -Q.71\d xy ) +0.71\d zx ). 
This lattice distortion brings about new terms C , D 1 
and E l in the Hamiltonian (see Appendix |A|) which pro- 
mote a mixing of the different kinds of orbitals. Thus, 
the uniform components of the orbital increase and the 
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hopping integral between the occupied orbitals r,j be- 
comes finite. As a result, the ferromagnetic ordering is 
relatively unstable in comparison with the AF spin order- 
ing due to the term A 1 . (2) The large JT-distortion with 
the Eg symmetry qeQe favors the orbital ordered state 
of any linear combinations of the d xy and d 
for A\ and A2 sites, and those of the d xy and d yz 



orbitals 
or- 
bitals for Bi and B2 sites. On the contrary, gT 2 QT 2 lifts 



the degeneracy uniquely as IV'A 



U\d x 



\d zx )), 



IVM 2 > = -^(\d xy ) + \d zx )), \%l> Bl ) = -±(\d xy ) + \d yz )) 
and \ipB 2 ) — ~^(\dxy) — \d yz ))- The obtained orbital 
state in the F spin phase in Fig. |](b) is the mixed 
state of the two. Here, the hopping integrals Ty's be- 
come finite for all NN bonds. Thus, a region of the fer- 
romagnetic phase shrinks by introducing gQ. (3) The 
LS coupling fixes the direction of spins in the ferromag- 
netic phase along [100]. The large orbital degeneracy at 
Xls = is lifted and one of the obtained orbital state 
is given by \i^ Al ) = \^b 2 ) = ^(i\dy z ) + \d zx )) and 



i xy i 



Although this orbital ordered 



state is included in the solutions at Xls = shown in 
Eq, (|3(]), there is no energy gain for the LS coupling at 
sites A2 and B\. This is because the F spin state is not 
compatible with the staggered orbital ordered state in 
the large limit of Xls- 
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IV. RESONANT X-RAY SCATTERING 

RXS was first applied to observation of the orbital or- 
dering in pcrovskite manganites.Ej By tuning the inci- 
dent x-ray energy to the K-e&ge of the transition-metal 
ion, the atomic scattering factor becomes a tensor with 
respect to the polarization of x ray and iSjiiensitive dra- 
matically to the local electronic structure.E 2 ] In this sec- 
tion, the scattering cross section of RXS in titanates is 
formulated. This is applied to the recent experimental 
results in YTiC>3. Let us consider the scattering of x ray 
with momentum ki , energy uji and polarization Xi to kf, 
Wf and A/. The electronic states at the initial, final and 
intermediate states of the scattering are denoted by \i), 
|f) and |m) with energies e^, £/ and e m , respectively. 
We start with the conventiona.1 form for the differential 
scattering cross section of RXSEjS 

d 2 a _ u f 



dVLdio f u>. 



A-f £|S| a *( £/ + (37) 



l/> 



where 



^ i (f\j- k i • S k i \ i \m){m\j kf ■ e kf x f \i) 



in 



£i - Em - W/ 

(f\jk f ■ e kf x f \m}(m\j-ki • e ktXi \i) 



iT 



(38) 



FIG. 5. The azimuthal angle dependence of the resonant 
x-ray scattering intensity, (a) I(100) aw , (b) J(001) CT7 r, (c) 
7(011) CTCT and (d) 7(011)^. The bold curves show the fit- 
ted results by the wave functions (a,b,c) = (—0.71,0,0.71). 
The doted curves show the calculated results by the wave 
functions {a,b,c) = (1/V2, 0, i/V2). The filled squares indi- 
cate the experimental data in YTiC>3 normalized by the in- 
tensity at (022) obtained in Ref. The inset of (a) shows a 
schematic picture of experimental arrangement. 



with the cross section of the Thomson scattering A = 
(e 2 /mc 2 ) 2 . T denotes the damping of a core hole, e k \ 
is the polarization vector of x ray and j k is the current 
operator. This form is rewritten by using the correlation 
function of the polarizability aip a as shown in Ref. |34|: 

J—.= A — Yj P 0'c*'P0c,n0'a'!3a(u>,i?), (39) 



dVldio f to. 



afla'f)' 



where 



-ri) 



(40) 



with K = ki — kf, uj = u>i — LUf and Pp a = 
(^kfXf )(3(eki\i)a- Now we express the polarizability oper- 
ator by the orbital operators Oir-y by utilizing the group 
theoretical analyses, aipa associated with the orbital op- 
erators at site I is represented by 



G 



E 

r 7 



(41) 



with (r 7 ) = (Eu 7 Ev), {T x x,Txy,Txz) and 
{T2X, T2y, T2Z). /r's are coupling constants which are not 
determined by the group theoretical analyses and Mr 7 's 
are matrices with respect to the polarization of x ray. Ex- 
plicit forms of Mr 7 are given by the Gell-Mann matrices 
A; introduced in Sec. || as Mr 7 = A; for (T7; I) = (Eu; 8), 
(Ev,3), Mr 7 - -A; for (T 2 x;6), (T 2 y;4), (T 2 z; 1), and 
M Tl = -i\ t for (T lX ;7), (T lV ;5), {T x z;2). As a result, 
the cross section for the static scattering in the orbital 
ordered state is obtained as 



with 



B(K) = Y j It 1 S Ti {0 Ti {K)). 
r 7 



(42) 



(43) 



Here, Sr 7 is the polarization part of the cross section 
defined by 



Sr 7 - ZkfXfMT-iSkiXii 

and (Or*f(K)) is the orbital order parameter 
1 



(°r 7 W) = ^E 



AK-fi 



(O 



;r 7/ 



(44) 



(45) 



with a number of the unit cell N. Equation (^2|) with 
Eq. (^3j) is derived by only considering the symmetries 
of crystal and orbital. Therefore, we do not specify, in 
this paper, the microscopic origin of the orbital depen- 
dence of the polarizability. However, it is reported re- 
cently that the mechanism of RXS based on the Coulomb 
interactions between 3d and Ap electrons well explains rel- 
ative scattering intensities at different reflection points in 
YTi03, rather than the mechanism based on the lattice 
distortions .c3E3 

The above results are applied to RXS in YTiC>3 where 
the detailed experiments have been done recently and re- 
ported in Ref. [2^. In the RXS experiments, the azimuthal 
angle scan, which is the rotational scan about the scat- 
tering vector K^-is crucially important to identify the 
orbital ordering]^ This is because the local symmetry 
of the orbital structure directly reflects on the azimuthal 
angle dependence. A schematic picture of the experimen- 
tal arrangement is presented in the inset of Fig. | (a). In 
this arrangement, Sr 7 in Eq. is replaced by 



S 



'fc/ A/ U{<p) VA/r-y V - 1 U~ 1 (y) eki \i j 



(46) 



where the matrix V describes the transformation from 
the crystallographic coordinate to the laboratory coor- 
dinate and the matrix U (if) represents the azimuthal 
rotation with the rotation angle ip. The GdFe03-type 



lattice distortion is also taken into account. A general 
formalism for the azimuthal angle dependence of RXS 
intensity is presented in Ref. Ejfjl The RXS experiments 
have been carried out at three orbital superlattice re- 
flection points of (100), (001) and (Oil), where we use 
the Pbnm orthorhombic notations. We fit the four sets 
of the experimental data: (100) with (A,, A/) = (cr,n), 
(001) with (er, tt), (011) with (a, a), and (011) with (a, tt), 
where a (tt) indicates the a (tt) polarization of x ray. 
From now on, the scattering intensity at (hkl) with po- 
larizations (Aj,A/) is denoted by I(hkl)\ i \ f . 7(100)™ 
and 7(001)™ are zero within the experimental errors. 
We assume that the orbital wave functions have a sym- 
metry of the Pbnm point group. Then, the coefficients 
Ci a in the orbital wave functions (see Eq. ((H)) satisfy 



the conditions Cah 



C 



B\zx 



Cf 



Caxzx — Ca 2 zx — Csxyz — Cb 2 vz = b, and Ca 1 xd — 
-Ca 2 x V = C BlX y = ~C B2 xy = c. In the case where 
a, b and c are real, the explicit forms of the scattering 
intensities are given by 

7(100)™ = 0, (47) 

1(100)™ = I Q (4I E E vSm<pcose) 2 , (48) 

7(001)^ = 0, (49) 

7(001)™ = 7 (47 T2 Ot 2 sinv?cos(9) 2 , (50) 

7(011)™ = I S(It 2 Ot 2 sin2<^) 2 , (51) 
7(011)™ = I 8(I T2 Ot 2 ) 2 

x (— cos2< / 9sin6' + sinip cos 9) 2 , (52) 

where the GaFe03-type distortion is neglected. 9 is the 
scattering angle, 7 = AN 2 /(2tt) and, Oe v and Ot 2 are 
the order parameters given by 



Ev = (0 Al Ev) = -^=(a 2 -b 2 ), 



(53) 



and 



Ot 2 = \((0 Ai t 2 x) - (0 A2 t 2V )) = ^c(a - b), (54) 

respectively. It is worth mentioning that 7(100)™ and 
other 7's arc reflected from different components of the 
orbital order parameters, i.e. Oe v and Ot 2 , respectively. 
From the above considerations, a, b and c are restricted 
so that c is finite and a b. We optimize values of a, 6, c 
and It 2 /Ie numerically within real numbers. It is found 
that the calculated results fit in the four sets of the ex- 
perimental data simultaneously when the wave function 
(a, b, c) satisfies the condition c > — a >> One of the 
best fitted results is obtained by a set of the parameters 
(a,b,c) = (—0.71,0,0.71), bijPig consistent with the re- 
cent experimental analyses J^jIIZI and It 2 /Ie = 0.45. The 
results are shown in Fig. |^ (bold lines) together with the 
experimental data (filled squares) . The fitting by the cal- 
culation is satisfactory. For comparison, we calculate the 
RXS intensity where the orbital wave functions are com- 
plex (dotted lines in Fig. §). We set a = 1/V2, b = 



7 



and c = i/y/2. In this case, the scattering intensities are 
explicitly given by 

7(100) CT7r = I (4I e Oev sm^cosd) 2 , (55) 
7(001) CT7r = Io(I Tl O Tl sin ^ cos #) 2 , (56) 

7(011) ffff = 7 i(7 Tl O Tl ) 2 

x (sin6> + sin^cos^) 2 , (57) 

and Tjo-'s are zero. Ot x is the order parameter for the 
magnetic dipolc moment defined by 

Tl = (0 AlTlX ) = ^ a * c ~ c * a )- ( 58 ) 

The discrepancies between the calculated results and 
the experimental data are remarkable, especially, in 
7(011) CTCr . This is because the order parameter Ot 2 van- 
ishes and Otx appears in this orbital ordered state. 
By using the orbital wave function (a, 6, c) 
0.71 . 0. 0.71 ^ obtained above, the spin wave dispersion 
relation is calculated by applying the Holstein-Primakoff 
transformation. The results are presented in Fig. |6| for 
several values of Rt- In the same sets of parameters, the 
Curie temperature Tq is calculated by the mean-field ap- 
proximation and calculated values are corrected by cop=. 
sidering the results in the high temperature expansion.EZI 
Rd = D zz /D xy , where D zz (D xy ) is the spin stiffness 
in the z direction (the xy plane), and Tq is obtained as 
(R D ,T c /J Tl ) = (0.75, 0.21), (1.00, 0.25), (1.23, 0.30) for 
R t =0.74, 0.84, 0.94, respectively. The spin wave in 
YTiC>3 is recently measured in Rcf. |3^ by neutron scat- 
tering experiments and its dispersion relation is found to 
be almost isotropic. A value of Jt x is estimated by fit- 
ting the experimental data by the calculated results of 
Rt = 0.84 where the spin stiffness is almost isotropic. 
The obtained value is Jt x = 10.9meV by which Tc is 
given by 32K. This result is consistent with Tc in YTiC-3 
of about 30K. 



V. SUMMARY 

In summary, we derive the effective Hamiltonian 
for spin and orbital states in perovskite titanates 
Ri_ x A x TiC>3. By taking into account the three-fold or- 
bital degeneracy, the orbital parts of the Hamiltonian 
are represented by the eight 3x3 matrices. The or- 
bital ordered states in the ferromagnetic phase for the end 
compounds have strong continuous degeneracy. This is 
owing to the three-fold degenerated ti g orbitals and the 
orthogonality of the electron hopping integral between 
NN ions, unlike the two-fold degenerated e g orbitals. In- 
troductions of the GdFeC>3-type distortion, the JT-type 
distortion and the LS coupling lift the orbital degener- 
acy and destabilize the ferromagnetic state. This implies 
that the high symmetry in the orbital space is favored 
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FIG. 6. The dispersion relation of the spin wave. Broken, 
bold and dotted curves are the results with 7?t=0.74, 0.84 
and 0.94, respectively. The cubic notation is used. The or- 
bital wave functions are chosen as (a, b, c) — (—0.71,0,0.71), 
and R.j is fixed to be 0.25. The ratio of the spin wave stiff- 
ness constant Rd = D zz /D xy and the Curie temperature 
T c is {R D ,T c /J Tl ) = (0.75, 0.21), (1.00, 0.25), (1.23, 0.30) for 
i?t=0.74, 0.84 and 0.94, respectively. 

for the ferromagnetic ordering. The scattering cross sec- 
tion of RXS is formulated by utilizing the same orbital 
operators adopted in the Hamiltonian. It is shown that 
the different components of the orbital order parameters 
are detected separately at the different reflection points 
and polarization configurations. The experimental data 
of the azimuthal angle dependent RXS intensities are well 
fitted by the present calculation. 

The present theory based on the derived effective 
Hamiltonian and new RXS formula is satisfactory to 
explain the several experiments about spin and orbital 
states in YTiC>3; the ferromagnetic structure, the az- 
imuthal angle dependence of the RXS intensity, the Curie 
temperature and the spin wave dispersion relation. It is 
thought that the spin and orbital states are controlled by 
the interactions between NN Ti sites through the virtual 
electron exchange processes and the associated JT-type 
lattice distortions. The present RXS studies suggest that 
the LS coupling may be irrelevant and the real orbital 
wave functions are realized by the inter-site interaction 
between orbitals. 

Beyond the static spin and orbital ordered states in 
i?Ti03, this Hamiltonian is applicable to the orbital ex- 
citations termed orbital wavenJ and to the hole doped 
systems Ri- x A x TiC>3. For the actual calculations, the 
present Hamiltonian represented by the 3x3 matrices 
Oir-y is more convenient than that by the 2x2 matrices 
Tu (see Eq. (p2h) with the constraint and the Hubbard 
type Hamiltonian. The study of the doped titanates and 
the orbital dynamics will be presented in a separate pub- 
lication. 
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E l = -V2(--\l-O l iEu )O jT2l 



V20 



2 12 

2 12 



,7j/ 1 3 



jEu 



(A7) 



By utilizing the pseudo-spin operator introduced in 
Eq. (122), these terms are rewritten as 



D l = 2(n ici T^+Tj x n jci ), 



(A8) 



APPENDIX A: EFFECTIVE HAMILTONIAN 
WITH A GDFEO3-TYPE LATTICE DISTORTION 



E l = 2{(n iai + n lbl )T l jx + T\ x (n jai + n jbl ) } . (A9) 



Effects of the GdFe03-type lattice distortion are in- 



cluded in the modified J term Tij introduced in Eq. ( |34[ ) 
The explicit form of this term is given by 



(Al) 



with 



{ij) V 



x i siR t 2D l + Rf [ ^-A l - C l + D l 



(A2) 



Ht 2 - -Jt 2 J2\J 



{ s, Hr2{D l + E l ) + R 2 t ( I A 1 +C' l + D 1 )} . 1 A:-! 1 



{ij) 



He = —Je \njrij - Si ■ Sj 



{ij) 

2 
; 3 



(A4) 



Ha! = — J Ax ^2 ( ~l ni7l 3 ~ &i ■ Sj 

{ij) V 



(A5) 



The parameter Rt = tc/t^ is caused by the GdFe03 



type distortion. AK B l , C , C 1 and D l are defined in 



Eqs. M), @, ©, © and dH), respectively. D l 
and i?describe new exchange processes induced by the 
distortion and are given by 



^o lT jl + \llo l 



3 V 3 



jEu 



(A6) 
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